Theory of the Diamagnetism Above the Critical Temperature for Cuprates 



O 

o 

(N 

< 

OS 



S3 
O 
O 

Oh' 



I 

C 

o 
o 



(N 
> 

in 

m 

in 
o 
m 
O 
-i— > 



i 

C 

o 
o 



X 



J.L. Gonzalez and E.V.L. de Mello 

Departamento de Fisica, Universidade Federal Fluminense, Niteroi, RJ 24210-340, Brazil 

(Received February 2, 2008) 

Recently experiments on high critical temperature superconductors have shown that the doping 
levels and the superconducting gap are usually not uniform properties but strongly dependent on 
their positions inside a given sample. Local superconducting regions may develop at the pseudogap 
temperature (T*) and upon cooling, grow continuously. As one of the consequences a large dia- 
magnetic signal above the superconducting critical temperature (T c ) has been measured by different 
groups. Here we construct a general theory to a disordered superconductor using a critical-state 
model for the magnetic response to the local superconducting domains between T* and T c and show 
that the resulting diamagnetic signal is in agreement with the experimental results. 
pacs74.20.-z,74.25-q,74.25.Ha,74.72-h 



I. INTRODUCTION 

There is an increasing interest on the effects of the mi- 
croscopic intrinsic inhomogeneities in high critical tem- 
perature superconductors (HTSC). The fact that these 
materials have doping level and superconducting gap 
that, even in some of the best single crystals, vary lo- 
cally inside a given sample, has profound consequences 
and a number of unconventional related phenomena are 
observed 1 specially in the normal region of the HTSC 
phase diagram. In particular, recent magnetic imag- 
ing through a scanning superconducting quantum device 
(SQUID) microscopy has displayed a static local precur- 
sor of the Meissner state at temperatures as large as three 
times the T c of an underdoped LSCO film 2 . Following 
up SQUID magnetization measurements on powder ori- 
ented YBCO and LSCO single crystals 3 ' 4 have shown a 
rather high magnetic response which, due to its large 
signal and structure, cannot be attributed solely to the 
Ginzburg-Landau (GL) theory of fluctuating supercon- 
ducting magnetization 5,6 . This is an unconventional be- 
havior because low temperature superconductors exhibit 
only a much smaller diamagnetic signal above T c which is 
believe to be a consequence of the thermal fluctuations 7 . 
Therefore such strong magnetic response was interpreted 
as due to the fluctuating diamagnetism produced by su- 
perconducting islands nucleated above T c 3,4,8 . On the 
other hand, the existence of superconducting islands 
above T c is a consequence of a non-uniform or disordered 
intrinsic carrier density, as pointed out by Ovchinnikov 
et al 9 . 

The tendency towards the formation of magnetic do- 
main lines or stripes was predicted long ago 10 and verified 
experimentally by different groups 11,12 . The stripes are 
formed by an instability towards charge phase separa- 
tion which means hole-poor regions of antiferromagnet 
(AF) insulating separated by hole-rich metallic domains. 
Recently, the microscopic intrinsic inhomogeneities in 
the charge distribution, consistent with the presence of 
charge domains or charge stripes, has been revealed by 
neutron diffraction 13 . A peak broadening in the atomic 



pair distribution function was measured and explained 
by a local microscopic coexistence of doped and undoped 
material 13 . Evidences for microscopic phase separation 
into normal and superconducting regions was also ob- 
tained by muon spin relaxation rate 14 . Therefore the 
inclusion of the intrinsic charge inhomogeneities or the 
tendency of the formation of charge domains in any cal- 
culations for the cuprates properties may be far more im- 
portant than previously anticipated. The main difficulty 
to incorporate such effect in quantitative calculations is 
that the exact form of the doping or hole distribution 
inside a given sample is not exactly known although it 
is a consensus that the local doping level becomes more 
uniformly distributed in the materials on crossing from 
the underdoped into the overdoped region of the phase 
diagram 1,13 . On the other hand, a considerable improve- 
ment in the understanding of the doping level and local 
superconducting gap functional form has been made by 
scanning tunneling microscopy/spectroscopy (STM/S) 
experiments 2,15-17 . 

Based on these findings, we have recently proposed a 
two-phase model (undoped and doped regimes present 
in a given compound) described by a bi-modal distribu- 
tion of the holes inside a given compound in order to 
model the charge distribution for a given compound be- 
longing to a HTSC family 18 . This distribution resembles 
a typical distribution of a spinodal decomposition of a 
binary alloy 19 . A two-phase model which arises from the 
fluctuation phase of the superconducting order param- 
eter was considered previously by Emery et al 20 . The 
basic ideas of our approach are: The undoped or hole- 
poor part of the distribution represents the AF domains 
and the hole-rich the metallic regions. The width of the 
metallic distribution decreases with the sample's aver- 
age doping since, as mentioned above, the hole distri- 
bution inside the compounds of a given family become 
more homogeneous as the average doping level or av- 
erage charge density increases 13 . Due to the spatially 
varying local charge density, it is also expected a distri- 
bution of the local T c 's, instead of a single and unique 
value as in usual metallic superconductors. Therefore a 



1 



given HTSC compound with an average charge density 
n m possess distributions of the charge density n(r), the 
zero temperature superconducting gap A (r) and the su- 
perconducting critical temperature T c (r), where the sym- 
bol (r) means a point or small region inside the sample. 
In this scenario we identify the largest T c (r) with the the 
pseudogap temperature T* of a given compound 21 . Dop- 
ing regions with n(r) > 0.05 have a metallic behavior 
with a decreasing T c (r) following a mean field pattern 22 . 
Upon cooling below T* the superconducting regions de- 
velop at isolated regions (mostly with n(r) < 0.05) as 
droplets of rain in the air and, eventually they percolate 
at the sample superconducting critical temperature T c at 
which superconducting long range order is established. 
This superconducting percolation scenario for HTSC has 
also been discussed previously but with somewhat dif- 
ferent approaches, by several authors 23,24 . Therefore, 
only at or below the percolation temperature T c a dis- 
sipationless macroscopic electrical or hole current may 
flow through the sample. Based on these ideas, we have 
used a bi-modal charge distribution in connection with a 
mean field calculation to reproduce the phase diagram of 
T* and T c as function of the average doping level for the 
BiiSriCaCuiO^+x and La,2- x Sr x CuOi families 18 ' 25 . In 
these chemical formulas x is the doping level which is 
similar to the average charge density, i.e., x — n m for the 
last (La series) and x = 2n m for the former (Bi) family 
of compounds 

The existence of the superconducting regions with dif- 
ferent T c 's in the material must manifest itself through a 
number of observable properties: we can number a few 
like the pseudogap phenomena, the resistivity and Hall 
coefficient temperature dependence, and so on. Here we 
want to discuss the normal state magnetization, whose 
study is the purpose of this paper. Thus, we develop 
a general theory of the magnetization of a inhomogc- 
neous superconductor which can be applied to the HTSC, 
in order to explain the recent anomalous magnetization 
measurements 3 ' 4 . In a similar fashion, L. Romano 8 in- 
troduced a theory for the magnetic response based on 
the fluctuations in the local order parameter which leads 
to a Kostcrlitz-Thoulcss transition and the fluctuation- 
induced diamagnetism 6 of superconducting islands above 
T c which reproduced well the details of the experimental 
data. Assuming that the superconducting islands started 
to be formed at T*, for measurements at T < T* and T 
close but above T c , some of these islands must be in a 
superconducting state and therefore they must yield a 
local superconducting diamagnetic response. We show 
below that the measured diamagnetic signal can be un- 
derstood within this picture of a superconductor formed 
by static domains with spatial varying T c (r)s. Since the 
local T c (r)'s can be much higher than the sample's su- 
perconducting T c , we apply the well known critical-state 
model (CSM) 7 to the magnetization response under an 
applied field. This approach is valid with the hypothe- 
sis of the measured temperature to be below that of the 
irreversibility line of the local superconducting regions . 



Some kink of hysteretic behavior was observed in YBCO 
and in LSCO samples 3 ' indicating that the irreversible 
temperature is above T c . We demonstrate here that this 
simple procedure is able to explain and reproduce the 
main features of the precursor diamagnetism measured 
behavior. 



II. THE MODEL 

From the above discussion, since the hole distribu- 
tion is not uniform, we can think of a HTSC sample as 
formed by different metallic an insulating regions of a few 
nanometers of magnitude. As a consequence of such dis- 
tribution, at a given temperature above T* there are in- 
sulating and metallic regions, below T* , there are insulat- 
ing, metallic and superconducting regions which started 
to condensate in the metallic regions at T* and below 
T c the long range superconducting order is achieved by 
a percolation transition while there still are some metal- 
lic and insulating regions present in the sample. As in 
Ref. 18 , we model these insulator and metallic regions 
by a bi-modal Gamma-type distribution with parame- 
ters which are based on the STM/S analysis 16 and which 
was used to derive the HTSC phase diagrams. 

Now we want to use the above ideas to discuss the 
magnetic response M{B) of such inhomogeneous super- 
conductors to an external magnetic field, between T c 
and T* . In this range of temperature, it is clear that 
only the superconducting regions or droplets with their 
T c (r)'s bigger than T c contribute to the diamagnetic sam- 
ple's magnetization. These superconducting region con- 
tributes to the sample's magnetization with a diamag- 
netic signal which depends on B and the overall sample's 
magnetization will be sum of all these contributions. 

In order to estimate the M{B) we follow the ideas 
and the procedures of the CSM to each superconduct- 
ing droplet. Upon applying an external magnetic field, a 
critical current (J c ) is established which opposes the field 
as J C {B) = a(T)/B according to Ohmer et al 26 . a(T) is 
a temperature-dependent local constant which we have 
taken to be proportional to (T c (n(r)) — T). Hereafter we 
call the hole local density n(r) of a given domain of n and 
the whole sample's average density of n m . For simplic- 
ity we take these superconducting droplets as cylinders 
of radius i?, which is sufficient small in order to have a 
constant charge density n and consequently the critical 
temperature T c (n) is the same within such cylinder region 
(As the temperature decreases, more droplets appear and 
the superconducting regions increases by aggregation of 
droplets of different n) . The CSM approach leads to the 
magnetic field dependence of the magnetization in each 
small cylinder as 26 : 

M 1 (B) = -—forB<B cl (1) 
Mo 

R 4R 3 o d5 

M 2 (B) = + — _ - ° ,B cl <B<B* (2) 

5hqB* 2 15fi B* 4 
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M 3 (B) = 
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(3) 



In the last expression, B* is the value of the applied ex- 
ternal field which produce full penetration inside a cylin- 
drical super conducting droplet, and it depends on its size 
as B* = y // 2a(T)noR- n o is the vacumn permissivity. 
The dependence of the local critical temperatures T c {n) 
is taken to be that of the pseudogap temperature of a 
given compound with its average charge density namely 
T*(n m ). Thus, T c (n) is a function which has its maxi- 
mum at n = 0.05 and decreases monotonously, that is, 
we take the local superconducting critical temperatures 
as a linear function of n, T c (n) = To — b * (n — n c ). To 
is the highest pseudogap temperature near the onset of 
superconductivity, i.e., n c — 0.05, and b is chosen in or- 
der to T c (0.3) = 0. Notice that, since T c (n) is a linear 
decreasing function of n, the droplets with lower values 
of n (just above n c ) are more robust to an applied field 
in the sense that its superconducting state is not easily 
destroyed by the external field. 

Since T c (n) is constant inside a superconducting cylin- 
drical droplet, the critical fields (B c i and B c2 ) inside the 
droplets will have their temperature dependence given by 
the GL theory, e. g., B cl (T) = B cl (0)[(l -T/T c (n)] and 
B c2 (T) = B c2 (0)[(l - T/T c (n)]. Taking into account the 
dependence of T c (n) on n, we arrive at the expressions 
for the critical fields B c i(T,n) and B c2 (T, n). A simi- 
lar functional form is supposed for B* due to the a(T) 
temperature dependence. Thus, 



B c i(T,n) = B c i(0)[l 



B c2 (T,n) = 5 c2 (0)[l 



B*(T,n) = B*(0)[1 - 
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When a given sample is submitted to an applied ex- 
ternal magnetic field B, the superconducting droplets 
with carrier concentration n for which the applied field 
is higher than their second critical field B c2 (T, n) — 
B C 2(0)[(1 - T/(T - b * (n - n c )], do not contribute 
to the sample magnetization. This condition is veri- 
fied for droplets with n > n max , where n max (B c2 ) = 
n c + T0/b- (T/b)/[l - B/B c2 (0)} is obtained inverting 
Eq.5. Since T c (n) is a decreasing function of n, only the 
droplets with n bigger than n max do not contribute to the 
sample's magnetization because their superconductivity 
is destroyed by the field B. Thus we expect that 



M(T, B) 



c {Be 



Where P(n) is the distribution function for the lo- 
cal hole doping level at the many clusters inside a given 
HTCS inferred in Ref. 18 . However, in the context of the 
CSM, depending on the intensity of the applied field, 
there are different possibilities in which each domain 
contributes to M(B). In the low field regime the su- 
perconducting clusters will contribute to the magneti- 
zation of the sample in three forms: there are some 
clusters, which are not penetrated by the field B, that 
is B < B c i(T,ri) and they contribute to the magneti- 
zation with perfect diamagnctism (Eq.l). These clus- 
ters have their carrier concentration in the interval, n < 
nc + T0/b- (T/b)/[l - B/B cl (0)]. The second group of 
clusters have their B c \ (T, n) lower than the applied field 
but B is also lower than B*(n, T). This group is partially 
penetrated by the field and they contribute to M{B) ac- 
cording Eq.2. These domains have their carrier concen- 
tration in the interval n c + T /b- (T/b)/[l- B/B cl (0)] < 
n<n c + To lb - (T/b)/[l - B/B*(Q)]. Lastly, there are 
some superconducting granules for which the applied field 
is higher than B*(T,n) but also lower than B c2 (T,n). 
These domains contribute to the magnetization accord- 
ing Eq.3. Therefore, for a sufficient low applied field, the 
general expression for the M(B) is given by: 



M (T, B) 



.(Bel) 



P(n)M 1 (B,T,n)dn 



rn max (B*) 

+ / P(n)M 2 (B,T,n)dn 

Jn ma x(B c i) 

/•nmax (B a 2) 

+ / P(n)M 3 (B,T,n)dn 



(8) 



Upon increasing the applied field the number of 
droplets that are not penetrated by the field decreases 
and for fields higher than B c i(n m ,T) there are not any 
granules contributing with perfect diamagnetism. No- 
tice that Bcl(n m ,T) is the maximum first critical field 
that a superconducting region can achieve. For B > 
B c i(n m ,T), all the droplets are partially or totally pen- 
etrated by the field and they contribute to the sample's 
magnetization according Eqs.1-3. 

rn max (B*) 

M(T,B)= P(n)M 2 (B,T,n)dn 

J nm 

rn max (B c2 ) 

+ / P{n)M 3 (B,T,n)dn (9) 

Finally for fields higher than B*(0)[(1 - T/(T0 - b(n m ~ 
n c )], all the superconducting granules of the sample are 
totally penetrated by the field and their contribution to 
the magnetization is 



P(n)M(B,T,n)dn. (7) 



rn m ax (B c2 ) 

M(T,B)= P(n)M 3 (B,T,n)dn. 

J nm 



(10) 
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This last three equations express the total M(B) for any 
value of the applied field. 

In order to obtain a reliable value of M(B) and to com- 
pare with the experimental results, we have incorporated 
the fluctuation magnetizations induced by the supercon- 
ducting order parameter, an effect which should be al- 
ways present, regardless whether the superconductor is 
more or less inhomogeneous. As noted in reference 27 , for 
superconducting droplets with a homogeneous order pa- 
rameter and with dimensions d approximately equal to 
the coherent length £(T), the Ginzburg-Landau model 
provides an exact solution for M fi U ct(B) . Here we use a 
simplified "zero-dimensional" for superconducting clus- 
ters of radius d smaller or near the coherence length ^(T), 
namely 3 : 



M fluct {T, B) 



2/5k B (ir£Td) 2 B 



4> 2 (T/T C - 1) + (ir£Bd)y5 



(11) 



where fcs is the Boltzmann constant, <I>o is the quan- 
tum flux and d w £ oc (T — T c )/T c . This last expression 
yields a linear Mfi uct (B) dependence for low fields and 
it has been incorporated in our calculations. The specific 
results for Mfi uct are shown in the insets of Figs. 1-2. 

Therefore Eqs.7-11 furnish the complete M{B,T) 
curve for a HTSC at temperatures which T* > T > T c . 
Our numerical results will be compared with the experi- 
mental results in the following section. 



III. RESULTS AND DISCUSSION 



The probability P(n) is also not well known but some of 
its feature was inferred from STM measurements and the 
phase diagram, as discussed in Ref. 18 . 



In Fig.l we plot the results of our model with the 
parameter which corresponds to a n m = 0.2. Accord- 
ing to our model this concentration corresponds to an 
LSCO ovcrdoped sample with a critical temperature of 
T c = 26.86K. In the calculations we have used B c i(0) = 
0.01T, B c2 (0) = 20T and B*(0) = 0.025T. The value of 
B C 2 is derived from experimental values 29 . The values of 
B c i(0) and B*(0) are not found in the literature, how- 
ever one can make a fair estimation of the ratio B c2 /B c i 
through the GL parameter n which is about 100 for some 
YBCO and LSCO compounds 30 . The inset shows that 
for T = 27. 8K there is a low contribution at low fields 
from the superconducting island in the sample while at 
high fields only a linear contribution comes from the dia- 
magnetic fluctuation signal. At T — 28K there is no more 
contribution from the superconducting droplets and all 
the diamagnctic signal is only due to the magnetization 
fluctuations. Note that the T = 28 K is the most week 
diamagnetic signal which appears in the Fig.l. The inset 
demonstrate that, at this temperatures or above, only 
the fluctuating magnetization is important. 



The above theory was developed to model the mea- 
sured magnetization curves of the Lai- x Sr x CuO<t family 
of compounds. The major difficulty is that the measure- 
ments of the pseudogap temperature curve for any fam- 
ily seems to vary depending on the technique used 21 . In 
fact there are differences even in the T c measured values 
for the same HTSC compound when reported by differ- 
ent groups, but such differences are small compared with 
those for the T* measurements. 

In the calculations we used values taken from the local 
precursor diamagnetic signal 2 and Nernst effect 28 mea- 
surements which display some local precursor order near 
T « 100K. This value is considerably lower than the 
values of T* from transport measurements 21 . A possi- 
ble explanation for this discrepancy is that the transport 
measurements detect the onset of pair formation while for 
the magnetic response a large density of pairs must be 
present in order to display a measurable local response. 
Thus, we have taken T = 80K, b = 320K for the carrier 
concentration linear dependence of the droplet 's critical 
temperature T c (n). In fact T c (n) is not necessarily linear 
but this is an approximation which fits reasonably well 
the experimental data but which may be improved in the 
future as new experimental data becomes available. The 
antifcrromagnetic cutoff for the carrier concentration is 
n c = 0.05 and the superconducting cutoff is n m = 0.3. 



The qualitative features of the measurements are en- 
tirely reproduced and are simply explained by our model; 
at low fields the perfect diamagnetism is expected for 
droplets for which the fields are lower than their B c \. 
We expect B cl to be weak because the superconducting 
regions formed above T c are small and isolated. By the 
same token, the droplets penetrating field B* should not 
be very strong what decreases rapidly the overall diamag- 
netic signal for field much weaker than B c i- As the ap- 
plied field increases, the magnetic response dies off and 
is reduced to the fluctuations. This is the reason why 
M{B) has a minimum at very low applied fields. In Fig.l 
we show this up-turn field near B 



up 



0.0001T which 

agrees with the experimental values 4 . It is worthwhile to 
mention that previously estimation for the up-turn field 
considering only the Lawrence-Doniach fluctuations 5 in 
a layered superconductor 3 yields expected values near 
B U p = 10T. These figures bring out the importance of 
the CSM applied to the superconducting islands above 
T c to explain the experimental results. 
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FIG. 1. Magnetization for parameters appropriated to the 
overdoped LSCO as calculated from Eqs. 7-11. The inset 
shows the change in the magnetization behavior as the anoma- 
lous contribution vanishes when the temperature is increased. 
This result is to be compared with the measurements from 
ref.4. 
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FIG. 2. Magnetization for parameters appropriated to the 
underdoped LSCO as calculated from Eqs. 7-11. The in- 
set shows the change in the magnetization behavior as the 
anomalous contribution vanishes when the temperature is in- 
creased. This result is to be compared with the measurements 
from ref.4. 



In order to study the effect of the local charge in- 
homogeneities, we have also performed similar calcu- 
lations with parameters appropriate to an underdoped 
compound. It is well known that, within the same fam- 
ily, overdoped samples have a more homogeneous charge 
distribution than the underdoped ones 13 . Thus, we have 
chosen a compound with almost the same critical tem- 
perature T c than that of the overdoped sample (shown 
in Fig.l) to single out the effect of the charge inhomo- 
geneities. In this simulation we used a carrier concen- 
tration of average concentration of 0.11 with its appro- 
priate distribution P(n) and with a critical temperature 
of 29. 9K. We have chosen the same field parameters, e.g. 
B c i(0) = 0.01T, B c2 (0) = 20T and B*(0) = 0.025T. The 
results from our simulations are shown in Fig. 2. In this 
case the overall diamagnetic signal persists for measure- 
ment temperatures higher than those for the overdoped 
sample. The B up field is dislocated to higher values but 
it is still much lower than the homogeneous Lawrence- 
Doniach 5 fluctuating field. This is a consequence of the 
fact that a underdoped sample has a more inhomoge- 
neous charge distribution than the overdoped compounds 
possessing regions where the field easily penetrates (with 
local n w 0.25) and others where are more robust to field 
penetration (with local n w 0.1 for instance). This fact is 
taken into account through the distribution parameters 
used in our simulation. The results for an underdoped 
compound are shown in Fig. 2 where can be also no- 
ticed in the inset that as the temperature approaches 
37K, the diamagnetic contribution of the superconduct- 
ing droplets decreases and for measurements temperature 
higher than approximately 37K there is only the fluctu- 
ations contribution. 



It is important to mention that, although our calcu- 
lations are very much in qualitative agreement with the 
magnetic imaging results 2 and with the powder oriented 
SQUID magnetic measurements 3 ' 4 , there are other exper- 
iments which have not detected a B up field 31 and have 
interpreted their results as solely due to the GL fluc- 
tuating diamagnetic theory 6 . However an inspection in 
their curves indicates that they have not looked in detail 
at the low field region. The absence of a B up field was 
also found on an optimally doped powder oriented YBCO 
compound 3 whose curves followed the pure fluctuating 
magnetic response. In this case, the possible interpreta- 
tion for a pure fluctuating diamagnetism signal, accord- 
ing to our theory, is that these samples must have a very 
good degree of charge homogeneity. This interpretation 
can be experimentally tested because these compounds 
must have T* almost equal or equal T Cl as may be the 

91 

case . 

Another point which hinders a good quantitative 
agreement with the experimental results is that the tech- 
nique used seems also to be very important in probing the 
details of the magnetization above T c ; while the scanning 
microscopy results on a LSCO film indicate the presence 
of a visible diamagnetic signal at temperatures as large as 
three times the sample's T c 2 , the SQUID magnetization 
measurements detected a diamagnetic signal on oriented 
powder only around 10% above T c 3 ' 4 . 

It is important to emphasize that, in our calculations of 
M(B,T), the islands that contribute to the diamagnetic 
signal are in the critical state, and the measured tem- 
perature T is below their local T c (r). In the other hand, 
the Romano's calculations yield an upturn field which 
is almost independent of the temperature while in Fig. 
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2 we can see clearly that it decreases with the applied 
temperature. Future experiments may distinguish which 
mechanism is more important or if some compounds ex- 
hibit one type and others exhibit the other. 



IV. CONCLUSION 

We have constructed a general theory of the diamag- 
netism applicable to a disordered superconductor with a 
distribution of regions with different local superconduct- 
ing T c 's. We have shown that the procedure reproduces 
the qualitative features of the unusual diamagnetic sig- 
nal above T c measured for several HTSC samples. This 
was done applying a CSM and the well known magnetic 
effects for fields between B c \ and B C 2 to the supercon- 
ducting regions formed at temperatures below T* of an 
inhomogeneous HTSC compound. The quantitative re- 
sults can be improved when experiments provide better 
estimates of T* and B c \. 

Our results demonstrated that, as concluded also from 
others different calculations 4 ' 8 , the measured normal 
state magnetization curves, the B up fields and the STM 
magnetic imaging results may be interpreted through the 
formation of static superconducting islands at tempera- 
tures above the sample's T c . 
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